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NOTATIONS 


The following symbols have been adopted for use in 
this dissertation: 


A 

a / a /a 

a-o 

d 

e 

®max 

®min 

F 

f 

g 

k 

L 

M 

ms 

mf 

N 

^s 

P,Q 

R 

RA 


= Viscous drag coefficient; 

= Position, velocity and acceleration co-ordinate; 

= Amplitude of soil vibration; 

= Diameter of the sphere; 

= Void ratio of the soil; 

= Void ratio at the loosest state of soil; 

- Minimum limit of void ratio of the soil; 

= Superimposed load on the sphere; 

= Frequency of soil vibration; 

= Specific gravity of soil; 

= Acceleration due to gravity; 

= Added mass coefficient of sphere; 

= A characteristic length; i 

= A non-dimensional constant, defined in Eg. (3. 6); 

= Mass of the sphere; | 

= Mass of the volume of soil displaced by the sphere; - 

= A non-dimensional constant, defined in Eg, (3,6) ; 

= A dim'-nsionloss constant equivalent to Stokes numberi 

= Non-dimensional constants, defined in Eq. (3.6); 

= Density ratio of sphere to soil mass; ; 

= Ratio of amplitude of soil vibration to the , 

characteristic length; 
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RB 

r 

S 

T 

t 


= Ratio of acceleration of gravity to soil acceleratic 
= Radius of the sphere; 

= A non-dime ns 1 one 1 constant, defined in Eq.(3.6); 

= Dimensionless time fector; 

= Time ; 


V 


Ratio of the non-dimens lonal velocity to the Stokes 
terminal velocity; 


= Stokes terminal velocity; 

V = Steady velocity of the sphere; 


w = Angular velocity of soil motion 

X,X,X = Dimensionless displacement, velocity and 
acceleration of thesphere, respectively; 


• (»♦ 
X, x,x 


= Position, velocity and acceleration co-ordinate 
of the sphere, respectively;; 


Xq = Half amplitude of the sphere motion; 

Z = Dimensionless velocity of the sphere. 


Greek Letters: 


B 

9s 

n 

^o 

fs 

^f 


= A coefficient to determine the effect of vibration; 
= Characteristic v lue of effect of vibration; 

- Unit v'/eijht of sohere; 

= Unit weight of soil 

= Ratio of acceleration of soil vibration to the 
acceletation due to gravity; 

= Threshol 1 of vibro-viscous state of soil; 

= Hass density of the sphere; 

“ ''ass density of soil 

= An interval, o 
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Coefficient of vibro-viscosity ; 

Kinematic viscosity of soil; 

Angle of internal friction of soil; 

Angle of internal friction of soil without 
vibration; 

Limit value of the angle of internal friction 
of soil duo to vibration; 

Stokes current function. 



SYNOPSIS 


The behaviour of Foundations, such as footings, piles 
etc. which are supported by che soil medium, when subjected 
to intense Vibrations either due to earthquakes (natural or 
artificial) , and other agencies, is an important aspect in 
Soil Mechanics. - ■4’he soil wh=n subjected to vibrations tends 
to liquefy and bthaves like a fluid as the intensity of 
vibration increases . The coefficient of viscosity ofthe 
fluidised soil medium (which is called vibro viscosity) is 
a significant property which influences the soil-structure 
interaction. Hence, to have a feel for tho behaviour of 
structures founded on such a soil environment analytical 
and experimental studies have been conducted in this 

y'' 

investigation 

In the study dascrib^-'d here in, the vibroviscosity 
of a loose-fine grained-cohosionless unsaturated soil 
deposit subjected to inuensi- vibration and its dependence 
on various soil parameters as well as the vibration 
characteristics has bcv^n studied experimentally and a mathe- 
matical relationship is established!^ 

^On liquefaction the soil behaves like a fluid and 1 

hone, the following Fluid Mechanics concepts have been 
brought h^-re to determine the soil characteristics as mentioned 
earlier. Stokes performed an approximate solution of the 
Navier-Stokes equations for a stationary fluid which several 
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investigators have attemDted to verify. Carstens solved 
the problem for an oscillating fluid making several 
simplifying assumptions. 

In accolorated motion of a submerged body, the fluid 
being accelerated exerts a force on the body. This force iSi 
represented by a product of the body acceleration and the 
fictitious mass called the virtual mass. Usually the 
Virtual mass is cither neglected or its magnitude is 

i 

determined from the irrotation-solutions for rectilinear 
motion, '^in the present analysis all the forces have been 
taken into account and solutions have been obtained using 
Ruhge-Kutta integration technique. ' 

-s^Iiab oratory tests have been conducted on the motion of 
a sphere in a soil medium subj'^cted to vibrations. The 
vibro viscosity ofthe soil has been obtained using the 
above mentioned analysis of the acccl^-rated motion of a sphe-^ 
in an oscillating fluid and has been compared with the values 
obtained from Stokes classical solution. The effect of grain 
size o:^thc soil, operating frequency, load superimposed on 
the sphere etc. have been studied. Results have been 
graphically presented and are compared with existing 
solutions . 



CHAPTER _ I 


INTRODUCTION 


1.1 General: 

Static investigations of non-cohesive soil show high 
resistance to deformation even if it iiS saturated and having 
low values of porosity. This is verified by field and 
laboratory tests i.e. load test, consolidometer tests etc. 

If any soil (cohesive or non-cohesive) is subjected to dynamic 
loads, the resistance to external loads decreases considerably. 
This IS verified by plunging vibrators, vibrating piles, pipes 
or cylinders with considerable diameter into the soil. 

The settlement of vibrating foundations or other 
structural elements are predominent in loose non-cohesive 
soils. Loose fine sands whether saturated or not show a loss 
of strength under the action of dynamic loads. The loads 
may come from various sources e.g. Earthquakes, Nuclear-blasts, 
machine foundations etc. The loss of strength results due to 
the changes in dissipative properties of soil e.g. forces of 
internal dry and viscous friction, forces of cohesion (which 
determine the initial resistance of soil to shear) , forces 
of external friction, hydrodynamic properties e.g, the 
coefficient of permeability and the pore water pressure and 
elastic and plastic characteristics such as Young's Modulus, the 
modulus of shear and the limits of elasticity and of plasticity. 
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It was also found, that due to vibration sandy fills 
and loose silty and clayey soils get compacted at first and 
if vibration intensity is increased more and more they loose 
their resistance to shear and a certain stage will reach when 
their mechanical properties are close to those of viscous 
liquids than of solids . The phenomenon is called vibrof loat- 
ation. The stage at which the soil behaves like a high viscous 
fluid is termed as liquefied stage and the phenomenon is 
called liquefaction. This phenomenon of liquefaction results 
in the failure of earth structures and foundations. This 
may include sliding of quay walls, settlement of earthen 
embankments ano filling up of mines with liquefied soil. 

1 .2 Imoortance and Purpose of the Study of Liquefaction 
Phenomenon: 

1.2.1 Studies in India; 

Importance of such studies is being felt in India 
because of the increase in construction activities in the 
seismically active zones. The Banas River Dam in Gujarat 
and Obera Dam in U.P. are two typical examples where the 
foundations are of loose sand and the dams are located in the 
seismically active zones. 

1.2.2 Importance from Field Evidences; 

During earthquakes, sand volcanoes and eruption of 

water and soil from the ground are often observed and especially 
23 

in the NIIGATA earthquake of June 16, 1964, extensive damage 
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to the structures occured as a result of decrease in bearing 

V 

capacity of soil. This decrease in bearing *apacity of soil 
IS due to two phenomena/ one of which is liquefaction of the 
soil which extends from the ground surface to some depth and 
the other is the loosening of the soil due to upward flow 
of water which resulted from liquefaction of the soil at 
some distance below the ground surface. 

Another typical example of liquefaction from earth- 
quake is TOKACHIOKI^^ earthquake of 1968. Due to this earth- 
quake saturated sands underwent complete liquefaction at 
NANAEHAMA beach, HAKODATE, HOKK/ilDO at coastal stretch of 
newly built hydraulic fill, resulting innumerous sand volcanoes 
of 2 ems to 50 ems high and the craters were 2 ems to 130 ems 
in diameter, grain sizes were ranging from 0.002 mm to 0.2 mm. 

It was observed at a plant in HACHINOHE city where a part of the 
building site was located on the back-fill sand prior to 
construction of the plant buildings, the loose sandi was 
densified by Vibrof loation and the main building was supported 
on the point bearing piles , the structure on piles and 
those on virgin ground showed little damage, where as the 
improved back-fill area underwent comolete liquefaction 

causing many structures to tilt settle or floatup. 

23 24 

In Ohsaki's ' opinion the following five items 
are the com itions under which the ground would liquefy during 
an earth quake of a considerable intensity, resulting in heavy 

y 
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damage to structures; 

1) Most of the soil from the ground surface down to a depth 
of 15 to 2o m consists of sands, 

2) It IS submerged under the ground water table, 

3) Its age of deposition is young , 

4) It consists of pure, uniform sands of medium gram size, 
i.e. i) the content of silt and/or clay is less than 10 
percent, ii) its 6o percent grain size is in the range of 
0.2 to 2,0 mm, and iii) its coefficient of uniformity is 
less than 5. 

5) Its comoactness is loose and N -values of the S.P.T. are 
smaller than the critical value represented by = 2 Z, 
in which Z denotes a depth in metres. 

During the earthquake of New Madrid in 1811, as quoted 
3 7 

by Karl Terzaghi , many slopes failed. Reason he gave was that 
the submerged masses of loose sand might under the impact of 
an eairthquake shock, temporarily assumed the character of 
suspension which flowed like a viscous liquid. 

1.3 Scope of this Investigation: 

The studies of Vibro floitation, particularly on 
liquefaction, were limited to saturated or partially saturated 
soils in English countries, although these tests were conducted 
for dry soils, in addition, by the Russians. 
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Linuef'ictxon studies/ which the present work deals a 
in particulnr, are being limited within the dry sands having 
Vory fine to medium grain size subjected to steady state tsfcertical 
Vibration. It has been tried to present a mathematical model 
for the sinking behaviour of asphere. It has been assumed 
that the soil flows like a viscous fluid with having very 
high viscous resistance, coefficient of which is known as 
vibro-viscosity . An expression is develooed for the acceler- 
ated motion of spherical bodies in an oscillatinj fluid, 
assuming the eiistance of a single degree-of -freedom system. 

Th© phenomenon of liquefaction results from the 
increase in pore water pressure and consequent reduction in 
shear strength, as it is being noted by several Investigators. 
Liquefaction studies in dry fine grained soils is important 
due to the fact that the surface of the ground is more 
vulnerable to liquefaction compared .fco the depth below. 

In finegrained soils the reduction in shear strength may be 
upto cent percent and as this strength reduces the soil behaves 
more anl more like a viscousf luid. 

The phenomenon of liquefaction of the fine sands 
have been studied both in the laboratory and in the field by 
various Investigators for moist soil. The factors that affect 
the liquefaction are: 

i) Grain Size Characteristics 
ii) Density of deposit 
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111) Surcharge weight 

iv) Vibratioa characteristics 

v) Location of drainage 

vi) Dimension of deposit etc. 

The vibration characteristics used here will be of 
steady state motion of simple Harmonic type, the other type, 
usually used, may be repeated loading. The soil property, at 
Its liquefied stage, is being characterized by the coefficient 
of vibro-viscosity , and the present work willbe to study the 
affect of frequency, amplitude and acceleration on the 
coefficient of vibro-viscosity , both theoretically and experi- 
mentally. This study will ultimately reveal the parameters 
which affect the viscous property of the soil continuum on 
vibrofloatation . 

Chapter II concentrates on the review of past literat- 
ure in the light of the present Investigation. 

Chapter III deals with the theoretical formulation 
of the subject matter and solutions. All theoretical 
derivations and results are presented in this Chapter. 

Chapter IV describes the experimental setup. A 
brief description of the each component and limitation of the 
ejq^erimental derivations etc. are reported here. 

Chapter V deals with the analysis of test results and 
its correlation with the theoretical findings and comments on 
the results performed. 
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The lest chapter (Chapter Vl) deals mainly with a 
general liscussionj. scone for further studies, commencs on 
the theoretical as well as exoerimental findings. 



CHAPTER II 


h GEMEHJ^JL, DISCUSSION AND LITER/.TURE REVIEW 

2.1 Introduction: 

PoTcrovsky (1934) and associates, as quoted by Barkan 
were the first to observe a decrease in strength c£ dry and 
moist sands under vibrations. They showed that the coefficient 
of internal friction depends on the Kinetic energy of vibration, 
as the energy increases, the coefficient decreases, approaching 
a value 25 to 3o percent smaller than that observed before 
vibration . 

In English known countries probably the first attempt 
to delineate conditions under which a considerable reduction 
in shear strength of cohesionless soil occurs, as quoted by 
Seed and Lee^^, was the 'critical void ratio’ approach suggested 
by Arthur Casagrande . It was noted that during shear, dense 
sands tend to expand where as loose sands tend to decrease in 
volume; thus, for any sand there will be some initial void 
ratio, termed the 'critical void ratio', for which no volume 
change during drained shear, (correspondingly, no pore-water 
pressure changes during undrained shear) , will occur. It was 
reasoned, therefore, that the sand deposits having a void 

ratio above the critical value and therefore tending to 

T 

contact during shear, would, under undrained conditions 
develop positive pore-water pressures that would possibly 
become large enough to produce liquefaction. Conversely, 
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deposits having an initial void ratio below the critical 
value would tend to dilate during shear, producing a decrease 
in pore-vjater pressures and a corresponding increase in 
effective stress, under undrained conditions, so that high 
strength and stability would be developed. 

2.2. Qualitative Approach; 

2.2.1 Tests with dry Soil: 

Barkan^ has carried out tests with dry sands and 
obtained a relationship between the acceleration of vibration 
and the coefficient of vibro-viscosity assuming the sand to 
behave like viscous fluid. He assumed the time versus sinking 
of a sphere to be linear and using Stokes law for static fluid 
he found out the coefficient of vibro-viscosity- Limitations 
of his studies will be discussed in the Article 2.4. 

2.2.2. Tests with Saturated Soil: 

Barkan^ also conducted tests with different moisture 
contents in sandy soils. He states, "The forces in Cohesion 
in soils depends on their moisture content which makes possible 
to assume that the coefficient of vibro-viscosity also depends 
on moisture content" . 

Maslov^^ gave the theory of * f liberation ' , according 
to which every sand can be subjected to consolidation under 
sufficiently intense vibrations. The intensity of vibration 
depends upon the state of compaction of sand. The compression 
of sand leads to considerable increase in pore-water pressure . 
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The sand grams may be looked upon as i£ they have lost their 
weight. Thus, depending upon the initial density of sand, there 
IS a 'critical - acceleration' value beyond which there will 
be an increase in pore-water pressure and decrease in strength 
and further consolidation of the soil. 

Seed and Lundgren^*^ could not liquefy a loose sand 
sample with transient loading, at a strain rate of 40 inches 
per seconJ, although an appreciable decrease in strength was 
observt-d in this strain range. They observed "Deformation 
of a donsv. sand causv^s an increase in volume, termed dilatancy, 
and deformation of a loose sand causes a decrease in volume. 

Thus if a loose saturated sand is rapidly deformed, the grains 
tend to move into a more compact condition, and for this to 
occur, water must flow from the voids. If there is no time 
for drainage to occur any load on the soil must be carried 
instantaneously by the pore water, hence, the pressure 
between the grams is considerably reduced, with a correspond- 
ing reduction in the shear strength of the soil. If the strength' 
IS reduced to zero, the sand is said to liquefy. On the other 
hand, when a dense saturated sand is deformed, the void 
ratio is increased, and if there is no time for seepage into 
the sand, the water can no longer fill the voids. As a 
result the pore-water is put into tension (some times called 
negative pore-water pressure) , there by increasing the 
pressure between grams and causing a temporary increase 
in strength" . 
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Matsuo and Ohara reported results of tests on loose 

saturated sands performed on a vibration table. They found 

a sudden increase in pore-water pressure occured at a definite 

acceleration. They confirmed the concept of 'critical 

19 

acceleration' given by Maslov. 

Florin and Ivanov^^ reported their vibration tests 
on a 20 ems thick sand deposit, both under steady state and 
transient vibration. They noted: 

i) If the sand is subjected to shock loading, 
the whole stratum get liquefied at the same 
time, while under steady state vibrations the 
liquefaction starts from the top and proceeds 
downwards . 

ii) The gram size and superimposed load affect the 

time for which sand remains in liquefied state. 

Ill) Presence of surcharge on the surface of the sands 

affected the mechanism of liquefaction. 

The above mentioned points were also raised by 
38,39 

Yoshimi 

13 

Huang-Wen-Xi obtained test data from triaxial tests 
performed by Shaking the sample on a vertically vibrating 
table. He used, density of sand sample, acceleration of 
Vibration and lateral pressure as variables in his study. 

He showed that the magnitude of pore-pressures depend not 
only on the physical properties of sand but also on the 
vibration characteristics. 
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The phenomenon of liquefaction/ he concludes, occurs due to 
loss of contact between the Grams by dynamic action, there 
by causing gradual transfer of shear stress from grain to 
pore-water . 

27 

Prakash and Mathur studied the effect of frequency 
and amplitude of vibration on samples of same cross-sectional 
area but of different heights. Their results show: 

i) The pore pressure attains a maximum value prior 
to resonance, 7i.t resonance the pore pressure gets 
partly dissipated. 

ii) Excess hydrodynamic pressure induced is larger 
than that is required to cause liquefaction. 

Ill) Liquefaction of top layers occur first and 
subsequently the bottom one. 

iv) Settlement of deposit varies characteristically 

with acceleration of table. 

3 1 

Seed and Lee concluded on the basis of their test 
results that the void ratio of the sand, the confining pressure 
acting on the sand, the magnitude of the cyclic stress or 
strain and the number of stress cycles to which the sand was 
subjected determined the danger of liquefaction of a saturated 
sand. Their conclusions were: 

i) Higher the void ratio more easily liquefaction 
will occur, 

ii) Lower the confining pressure more easily liquefactio 
will develop. 
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Hi) Larger the stress or strain, lower the number 
of cycles require to induce liquefaction. 


2 .3 Analysis of Vibro-f loatation? 

2.3.1 Effect of Acceleration on Soil Behaviour: 

Barkan^ made extensive studies and gave an empirical 
relationship between the coefficient of internal friction 
of sand and acceleration of vibration. It shows as acceler- 
ation increases the internal friction decreases and the 
decrease is rapid upto an acceleration ratio of 5.5 (Fig. 2, 6) 
and the relation is expressed as; 


where, 


tan 0 = ( tan 0g^ - tan 0oo ) Exp (- ) + tan 


tan 0^^ = Value of the coefficient of internal 

friction without vibration, 
tan 0OO = Limit value of the coefficient of internal 
friction. 

= Ratio of acceleration of vibration to 
the acceleration due to gravity, 
p = A coefficient that determines the effect 

of vibration (for dry medium-grained sand, 
p = 0.23) 


It IS to be noted that at liquefaction there should not be 
any internal friction theoretically, but frckn the above 
evidence it will not be zero in any case. 
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2.3.2. Effect of Frequency and Amplitude on Soil Behavior: 

These, too were given by Barkan^. Fig. 2.1 shows 

thfc relation between the coefficient of internal friction, 

0 , and the amplitude of vibration for dry medium grained 
sand frequency remaining the same. The graohs show that the 
coefficient of internal friction of sand decreases continuously 
as the amplitude increases and at very high amplitudes it 
becomes asymptotic. 

The dependence of tan 0 for same angular frequency 
IS more complicated. As frequency increases upto 180 
radians per second, the coefficient slowly decreases; then, 
as the frequency increases still (from 180 to 250 radians per 
second) , the coefficient of internal friction sharply 
decreases; and subsequent increase in the value of frequency 
ha(S almost no influence on the coefficient of internal 
friction. The trend is shown in Fig. 2.2. 

2.3.3. Effect of Moisture Concent and Gram Size: 

It was shown by Barkan^ that the coefficient of 
internal friction decreaseswith the moisture content, the 
lowest value occunng at 13 percent of moisture content. 

He conducted tests on four varieties of sand under 
two regimes of vibrations to study the effect of grain size. 

If the effect of vibration is characterized by, S'/ 

tan - tan 0 

tan 0gt 


where ^ ^ 
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the relationship is given as shown in Fig. 2.3. 

2.3.4. Effect of Acceleration on Vibro-viscosity : 

Barkan^ was the first to study the effect of 
acceleration on vibro-viscosity . He considered the behaviour 
of finely grained sand subjected to vibratory loads to behave 
like a fluid and observed the sinking of a sphere in the 
vertically oscillating dry sand. He utilized Stokes law, 
which establishes the dependence of velocity of the motion 
of a sphere m a static viscous fluid on the resistance acting 
thereon, the radius of sohere and the coefficient of vlscoslt^ 
of the fluid. He expressed the relationship for a sphere 
falling un ler the action of gravity in a stationiary viscous 
fluid as: 

= I ( '^s - 


where, ^ = Coefficient of viscosity, 

V = Steady velocity of the sphere sinking, 
r = Radius of the sphere, 

tnit weights of the sphere and the 
fluid respectively. 

In obtaining , Barkan did not consider the 
vibration of fluid (here soil) at all but later, based on 
e> 5 )eriments he gave a relation between ^ and accleration of 
Vibration as (Fig. 2.4) 
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^ = a (1 - 'Iq) 

where, a = A constant coeffxcient, 

=:R''.txo of Acceleration of soil to acceleration 
due to gravity. 

= The threshold of the vibro-viscous state of soil. 
It has b^en established experimentally that the acceleration of 
vibration lower than = 1.5g, docs not, practically, affect 
the coefficient of vibro-viscosity,. where, g is the acceleration 
due to gravity . 

2.3.5. Effect of Acceleration on Porosity of Sand: 

Barkan^ describes this relation-ship for fine to 
medium grained sand on experimental evidence (Fig. 2.6) as: 

® = °min + ‘®max - 1+ loO 

where, 

e =Void ratio at acceleration ratio, 

®max — ''^oid ratio at the loosest state of soil 
®min ” Minimum limit of void ratioof the soil 
A = A value depends on moisture content to a 
maximum value between 0.82 to 0.88. 

2.3,6 Effect of Surcharge and Confining pressure: 

Florin and Ivanov^^ stated that the saturated soil 
which was placed at a depth more than 15 metres below the 
ground surface was hardly liquefied. 
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1 

Ohsaki, as quoted by Kishida'^ , showed in his report 
on damage to R.C.C. Building in Niigata that the upper surface 
of the liquefied soil layer in the most severely damaged area 
IS mostly situated at depths less than 8 metre. 

Kishida'*’ estimated the effective overburden 
pressure of liquefied soil during Niiqata earthquake to be 
aporoximately 0.8 kg to 1.8 kg. per squared centimetre. 

He states that a sacurated sandy soil is not likely to liquefy 
if the value of effective overburden pressure exceeds 2,0 kg/cm . 
He co-relates the 'N' values of Standard Penetration Test with 
effictiVG overburden pressureand shows the probability of 
liquefaction of saturated sand as given in Fig. 2.7. 

2.4. General Discussion and Develomient of the Present 
Investigation: 

It can be noted that the critical void ratio is not 
constant for a given sand, but it depends on the confining 
pressure to which the sand is subjected... . This shows an 
inadequacy of the critical void ratio concept given by Arthur 
Casagrande, as brought out in section 2.1, 

Seed^*^ writes, "Because the dilation tendencies 
are smaller at high confining pressure, the critical void ratio 
decreases as confining pressure increases. Thus it has some 
times been concluded that a saturated sand at any given density 
IS potentially less stable under a high confining pressure 
(producing compression characteristics) than under a low confin- 
ing pressure (Producing dirtant characteristics) " . 
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The afforementioned approach, by Arthur Casagratide 

along with Seed's comment, can provide a valuable guide to 

the behaviour of saturated sand for which the critical void 

ratio is established, that is, single application of stresses 

building up from zero to some maximum value* Casagrande 
30 

noted, as Seed quotes, in presenting the concept of critical 
void ratio, volume changes under cyclic loading conditions 
are quite different from those occuring under one dimensional 
static loading condition, and it could hardly be ejqjected that 
the critical void ratio concept would be applicable to vibratory 
loading conditions. 

19 

However, Maslov , noted the ’’breakdown of 
Structures associated with sand where the porosity is below ’ 
'critical* and, on the contrary, the irreproachable behavior 
of many structur-s with the sand porosity above critical”. 

The afforementioned inadequacies lead the several 

27 19 

engineers like, Prakash and Mathur , Maslov etc. to 
attempt to establish the conditions producing liquefaction 
interms of acceleration at which the phenomenon can be observed 
to develop. 

The limitations of Barkan' s study as shown in 
sections 2,2,1 and 2,3.4, is due to the fact that he did not 
consider the vibration of the soil in determining the coeffi- 
cient of vibro-viscosity and at the same time he assumed 
sinking of the sphere with time to be linear, which is not 
true as it is evident from his e^qjerimental curves , 
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It xs, therefore, proposed to consider the vibrating 
soil as vibrating viscous fluid and then findout a solution 
for the sphere in an oscillating fluid under the action of 
gravity . 

2.5 Accelerated Motion of Sphere in a Fluid; 

Stokes investigated the simple harmonic and 
rectilinear oscillations of a sphere, a cylinder and an 
infinitely long flat plate in a viscous fluid. He omitted 
the convective acceleration terms in the Navier-Stokes 
equation and derived the expressions for forces exerted by 
the fluid on these objects. Each expression consists of two 
terms, one involving the acceleration and the other the velocity 

and both includes viscosity. 

2 3 17 

Later Basset ' , Boussinesq and Oseen studied 

the rectilinear motion of a sphere which has a rapid but 
arbitrary acceleration in a viscous fluid. They also 
omittcl convective acceleration terms in the Navier - Stokes 
equation in deriving their equation of motion. They postulated 
that the force on the sphere depends not only on its instant- 
aneous velocity and acceleration but also, on a term, which 
will be termed as Basset’s ’History Integral*, which represents 
the effect of its entire history of acceleration. Each effect 
is represented by a separate term. It is also important to 
note that the acceleration term does not include viscosity. 
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Infact/ this term is the same as the expressian derived for ar 

inviscid an-l irrotational flow. 

2/3 

Basset ' presented the equation of motion 
(accelerated) of a sphere falling in a viscous fluid under 
the influence of gravity in which the resistance is a linear 
function of velocity. He also assumed that for slow motions 
(as IS the case here) the squares and the products of 
velocities can be neglected. Thus the equation of Basset reads'. 



= (mg - m^) g ,,, (2.1) 

where m^ = Mass of the sphere 

m^ = Mass of the fluid displaced by sphere 
k = Added mass coefficient ( = 1/2 for sphere) 
v_ = Velocity of the sphere 
d * Diameter of the sphere 
^ = Density of the fluid 
^ / IS any time/ 

^ = Coefficient of viscosity of the fluid 
g = Acceleration due to gravity. 

The solution for equation 2.1, is given by severaJ 
persons like Basset^^^/ Lamb^"^/ Brush et al ^/ Hjelmfelt etial 
^^/ Lord Rayleigh^^ etc., but each solution has go^ 


Oder etal 
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limitatxon one orother and also they can not be applied to the 
case of an oscillating fluid. 

Based on these concepts, a detail out line of the 
solution for the accelerating motion of a sphere in an 
oscillating fluid will be presented in the next chapter. 



CHAPTER _ III 


3.1 I nt roduct ion ; 

It has been established m Chapter II, that the 
soil, subjected to intense vibrations or shock loads, tends 
to behave like a fluil. In order to study the behavior of 
structures on such a soil, some of the parameters, like the 
vibroviscosity of the soil etc. have to be analysed, With 
that in view, the classic problem of the accelerated motion 
of a sphere in vibrating soil continuum has been analytically 
studied in this chapter, idealizing the soil medium to 
behave like an oscillatingfluid. Numerical solutions have 
been attempted for the resulting equation of motion and 
^results are graphically shown. 

The basic study of the accelerated motion of a sphere 

0 

in an oscillating-fluid was due to Carstens who obtained 
the solutions making several simplifying assumptions. In this 
study the basic equations have been modified in the light of 
the analysis presented in sec, 2.5, to suit the problem 
under consideration. 

3.2 Accelerated motion of a Sphere in an Oscillatingfluid; 

It will be assumevd that the sphere le falling 
under the action of gravity. An one-degree of freedom 
system will be assumed throughout the analysis and it is 
also assumed here that the fluid is oscillating with 
Simple Harmonic Motion. 
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Now let, 

a = Position CO— or'-iinate of the fluid, 

a = Velocity co-ordinate of the fluid, 

a = Acceleration co-ordinate of the fluid; 

X = Position co-ordinate of the sphere, 

X = Veclocity co-ordinate of the sphere, 

X = Acceleration co-ordinate of the sphere; 

rOg = Mass of the sphere, 

C = Mass density of the sphere, 
s 

= Mass of the volume of fluid displaced by the 
sphere , 

= Mass density of the fluid, 

ytc r'f = Coefficient of vibro— viscosity and Kinematic 
Viscosity of the fluid, respectively, 
k = Added mass coefficient of the sphere 
= 1/2 for sphere* 
d = Diameter of the sphere. 

The forces involved in this type of motion are 
as follows: 

mg X = Inertial reaction of the sphere, 

mf a = Force exerted by the fluid on every fluid volume, 

k m^Ca-x) = Virtual mass effective force. 

Dissipative forces on the sphere o 

t V ( a-x) df 

o o i d?* 

= a-n-AdC^-x) + 4 dVTT f. \ f..:. 

'' 2 Jo 

This IS directly from Eq. 2.1. 
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On account of the pressure gradient in the fluid, 
if the fluid IS accelerated in the direction a, there is a 
force a) on the sphere. The force T<mf(a-x) is the virtual 

mass multiplied by the relative acceleration between the fluid 
and sohere. This simole treatment is possible in a system 
with one -degree-of -freedom, since the acceleration is a true 
vector in this co-ordinate system. 

The History Integral: 

This term is a type oi' histoiry term which indicates 
that the residtance to motion at any time t, due to unsteadiness 
of the flow, is, a apart, a function of the resistance at a 
previous time. This corrects for the transient character 
of the velocity distribution. The derivation of the appropriate 
History Integral for the oscillating fluid has been derived 
in Appendix -II. 

Thus using Newton' s law, the equation of accelerated 
motion of a sphere in an oscillating fluid can be written ms; 

mg X j= mf a + (a-x) + (mg-mf) g + 371^ d (a-x) 

,t (a-x) dr 

( dT^ ... (3.1) 

Jo V't -r 

Assuming the fluid motion as; 
a = a^ Sin wt 



. .. (3.2) 
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where , 

o 

w 


= Half araplit-ude of fluid motion^, and 
= Angular vf’ioc’iy of fluid motion. 


The Eq. 3,1 takes the fora; 


18 

X + 


X + 9 


+Kk)d 
^f 


2 es 


\/>V7r I 


t — (~ - w ao Cos wHd^' 

^ d^ dr 


(^s_ + 3^) ^ 

^f ^ 0 


A -r 


( i) 

2 (1+k) ^ 

- w a- Sin w t + -' 7 r= y 

-f 


18 '? 


+ ,1 


if - + k) d^ 
^f 


w Cos w t 


.. (3,2) 


Non-dimensional form; 


using; 


Eq, 3.3 can be written in the non-dime ns lonal form. 


X = 


X 

L 


T 


>?t 

d2 


^ X = XL 



where. 

Thus: 


L = Any characteristic length; Here the depth of 
fluid upto which the teadings can be taken. 

dt ^2 dT ' dt^ 
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Agairi/ let, -p- ~ R, the densit7 ratio of sphere to fluid, 

'■f 

and RA = ^ , RB = g/w^L 

Li 

= A number equivalent to Stokes number 


')? 


wd^ 


So then, the non-dime ns lonal form of the Eq. 3.3, reduces to; 


d^X 

dT^ 


18 1+ k RA 

(R+k) dT ^4 

o 

Cos (t/n„ 2) _ jR-rlL RB 


(R+k) N, 


N/^/TTd ? 

(R+k) 


%- ,4 

d 

X T7 


, , 2. 18 RA 

Sin (T/N ) - . 

s (R+k) 


— ( Cos w )d7> 

dr dr ° '■ 




(3.4) 


1/2 

Multiplying each side by (t-D and integrating with respect 
6, 29 


to ^ 
•t 


, denoting the left had side as F(T) ; 

^ d /dx _ a_Coswrj 


F(T) 


df = 


9sJ 


FVl 

/TTd^ 


(R+k) 




'o 


f 


dh dr 


> r ^ 


dr 


(3.5) 


r 


Let G(t) = 


- w a Cos w r ) df / ^ 


0 ^ r -d- 

0 




F' (r)dr 


-r 


^ s / t-r 
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Then the double inteqral, 


( G(t) d'?- ^ using ‘Abels' theorem; 

J ^ >1 t-T 


0 


= -rr [i(t) - F(o)l 


Thus Eq. 3.4 yields/ 


0 


F(T) 

4^ 


df 




TT d^ 


(R+k) 


^pL 'it 


Using (x - a) 


w a. 


= O 


t = 0 


Thus the non-dime ns lonal form of Eq. 3.1 reduces to; 

d^X , „ ^ ^ + Q Sin (T/Np) - s Cos (T/U, 

— r + P dT ^ pr; dT ' s 

dT^ ^ 


- M - pz: (Cos (T/Ng^) -l) = O 

■nTt 


Where, P = 18/(R+k) 


S = P. RA/N^ 


O - .RA/S M = f-T • — 

o 


N = .442 P. RA/N, 


Cos wt) 


(3.6) 
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The Eq. 3.6 is the final equation governing the 
accelerated motion of a sphere in an oscillating fluid in the 
dimensionless form. 


3 .3 Numerical Solution for the Governing Equation: 

A fourth order Runge-Kutta method will be used for 
solving the Eq. (3.6) . 

The equation is broken up into two first order 
simultaneous differential equations as follows. 


6X ^ 
dT 


, (3.7a) 


^ = S* Cos(T/N ^) +M + —r (T/Ng2)_i) 

dT >/ T 

- P (1+ .442/s/~T )* Z -Q* Sin (T/N ^) 

... (3.7b) 


Thus, there are two dependent varibles X and Z and 
one independent variable T. 

Difficulty \vith the Eq. (3 ,7) is, althouglo, 


but 


dT 


d^X 

2 

dT 


= RA/N^^ and X 

o 


T=0 


= 0 , 


T=0 


dZ 

dT 


, can not be evaluated directly 


T=0 


and only a limiting value is possible at T = 0, owing to 
Its singularity at T = o. 
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However, the limiting value of ^ as T — ^ 0 can be found 

as, — = M, vjhen T if very close to zero. 

dT 

Upto a certain range of T the value of X and Z will be 

Z 

obtained from the limiting value of and after that, 

the calculation will be carried out using Eg. (3.7b} . The 
validity of such an approximation is verified in the sec. 3.6. 

3 .4 Runae-Kutta Method with Gill's Variation: 

This method comes under the general category of 

Runge-Kutta methods and is a modified version of Runge-Kutta 

fourth order method. The derivation of this method can be 

2S 

found in many text boohs on numerical analysis 

The Runge-Kutta numerical integration method is a 
non-iterat ive, step by step and self starting procedure. This 
method is suitable here as the problem dealt here is an initial 

P 

value problem having two boundary conditions at the starting 
point of the integration. mAccuracy of this method is derived 
by using several estimates of the dependent variable for each 
increment of the independent variable. Usually the method is 
set up for first order equations . However, it can be easily 
adopted to higher order equations. A second order differential 
equation, as in this case, is handled by first transforming 
them into a set of first order differential equations (Eq.3.7). 

The truncation error for one integration step is of 
order (dT) ^ for this fourth order method, where (dT) is the 
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iticrement of independent variable. Theoretically/ there are 
no limits on the increment size in regard to convergence and 
stability. The size of the time increment *f integration 
is decided so that it is not too small as to require excessive 
computation labour or result in large round off error; neither 
It should be so large as to give rise to large truncation 
errors (the variation with respect to step size is shown 
in Fig . 3.2) 

When the time increment of integration is increased 
less computer time is needed for the solution of a given 
problem/ but beyond a certain limit of increment the solution 
obtained would deviate more and more from the true solution, 
and instability in the numerical integration is said to have 
occured and so there is a limit as to the larger time increment 
that can be used in the integration. The increment (dT) is 0-05 


as is taken here. 
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3 . 5 Conparision of the Solutions Obtained with the Existing 
Results ; 

8 

3.5.1 Cars teas ’s Solution: 

This solution IS essentially the same equation as 
3.1, neglecting the history integral of Basset,. The non- 
dimensional form of the equation is thus; 

j 2 v" o 0 

= S* Cos (T/N^ ) - P qr- - 0 Sin (T/N + M 

dT s 

... (3.8) 

Results of the equation (3.8) are shown in Tigs. Sitosir 
and discussed in the article 3.6. 

3.5.2 Wagenschein' s Solution: 

O 

As quoted by Carstens Wagenschein conducted some 

experiment on the oscillating motion of sphere m an 

oscillating fluid. Carstens checked these results using viscous 

17 

drag and Added mass coefficients given by Stokes, using Lamb's 
solution . 

The maximum displacemen'*" of the half amplitude of 
fluid motion to sphere motion in an oscillating fluid 
given by Konig (quoted by Carstens) is as follows; 


Xq W (1 + k) 2 + (A/mf w) 


... (3,9) 


■"o 


(R+ k) ^ + (A/mfw) 
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where/ = Half amplitude of sphere motion and fluid motion 

IS, a = a Sin wt, same as shown in equation (3.1' 

^ o 

To apply this to the case of oscillating sphere/ the 
coefficients A and k have to be changed to: 


Viscous drag coefficient/ A = ITT^dC 


k = h + 


... (3.10) 


Iv 

where, N = ^ 


For such a case the equation (3.5) reduces to. 


p* (1 . + ,442AA? ) * (( 8/n/) * Cos (STA^g^) - 


^i±3sl_ ^ (64/N ^) * Sin (8 T/N ) 

(R+k) s s 


.. (3.11) 


Where X = x/a^ , = half amplitude of fluid motion 


P = 18/ (R+k) 


N 




8 y 

wd^ 


These results are compared using, A - 1 .0 and 
k, the added mass coefficient to be equal to half. They 
were plotted m Fig. 3.14 and are discussed in the article 3.6 
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3 . 6 Discussion; 

To start with the solution of Eq. (3.6)^ the applica- 
bility of the Runge-Kutta integration technique to the 
present problem is justified by comparing the numerical solution 
of the static case of Eq. (3.1), with the existing solutions. 

IP 

One such solution is reported by Hjelmfelt et al and his 
solution IS verified using the following computations; 

Eq. (2.1) , can be rewritten as; 

TTfli ( + k ff) t. = ( P^) g - 3Tr/, ^ vd 

1 ^ (3.12) 

o 

which ultimately converts (Using Abel's Transformation) to; 


A ^ + (l+.443/>/ T) V = 1 


... (3.13) 


where, A = (R + k)/18, R = Density ratio 

V_^ = Stokes terminal velocity or terminal Reynold's 

number = (R-l) gdVlSf^ 

T = Time factor = t/d^ 

V = ) — / where v is the velocity of sphere, 

r Vj_ 


The solutions of the eq. (3.13) are shown in the 

Figs. 3.12 and 3.13, nhey show an excellent agreement with 

12 

the existing solution given by Hjelmfelt et al 
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The figures 3.1 to 3.11 show the solutions of the 
Eq. (3.6), depending upon the various magnitudes of non- 
dimensional variables namely, Ng, RA, RB and R. In most of 

O 

the cases Carstens's solutions are also presented and it has 
to be born in mind that the solutions reported by Carstens® 
are essentially the solution ofthe Eq. 3.6 neglecting Basset’s 
history integral. It is worth mentioninq that Carstens’s 
solution gives much higher values of displacement (as high as 
20%) than the actual solution (Eg. 3.6). 

To read the graphs the following notations are 
followed: 

In Fig. 3.1: A, denotes the proposed solution of the 
Eq. (3.6) , 

B, denotes Carstens's solution, and the 
density ratios are shown in the bracket. 

Here, B(5.0) denotes, the Carstens's solution with 
density ratio, R = 5.0, and same is the case for: A(5.0), which 
denotes, the proposed solution with R = 5.0. 

In almost all the results two variables were kept 
constant while the other two were varied. In Fig. 3.3, A(5.0)3 
denotes a true solution v;ith the density ratio, R = 5.0 and 
keeping R/. and RB constant, the magnitude of the other variable 
IS 0.1, which is the case 3 ox Ng variation. This procedure 
IS followed throughout the representation of the solution. 

As N_ increases, it can be observed that the deviation 
s 

of the Carstens's solution from the true solution is more. 



35 


The solutions are given in terms of displacement - time 
relationship/ keeping the practical utility of these solutions 
in view. 

In Fig. 3,3, it IS observed that the time- displacement 
relationship is a linear after a certain time interval, which 
indicates a little change in the nature of solutions with the 
change of kinematic viscosity. 

It Cc n be seen from the Fig. 2 .^ that as the amplitude 
of vibration vanes, the nature of solution also changes. As 
an example, for R/I = 0.01, alter a certain time the magnitude 
of velocity decreases an' then becomes steady, where as for 
the other valu's of Rk, the variation of velocity is of 
increasing order and then becomes steady. In lower density 
ratios, of course, the nature of variation is nearly same 
as for higher RA values. Thus, it can be concluded that at 
lower density ratios the role of amplitude variation is 
appreciable. This can be noticed also from Figs. 3.6 and 3.8. 

Excerpt for lower density ratios, the role of accelerat- 
ion ratio, R3, is not predominent. This can be observed in 
the Figs. 3.5, 3.7 and 3.9. 

It can be notice! from che Figures 3.10 and 3.11 
that at the higher valu-s of the variables (Ng, R/j., RB), the 
time - displacement relation is almost linear. 

The Fig, 3.14, gives a comparative study of the 
presented solution wich that of Konig's results. It is to be 
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noted that for higher density ratios, two solutions 
converge. Again, it is worth mentioning that at lower 
values of h , the variation is appreciable, again at higher 

o 

valuas of Ng(i.e. of V or lower values of w) the proposed 
solution deviates appreciably from Konig's solution. 

Althoug, Konig’s solutions are experimentally verified by 

3 R 

Carstens , Wagenschein / still, they are aooroximate 
solutions of Eg. 3.6 neglecting history term of Basset. 



CF PTER _ IV 


E XPERII'^ENT/iL EET U P 


4 . 1 Ob lective ; 

The present laboratory investigation was carried out 
with the follovJing objective^; 

To e'^^alrate the vibro-viscosity , yw. , using conventioKral 
m thod and its variatioi with respect to frequency, load and 
Grain Size che uac-erist i. :s of the soil, and then compared 
with the solution., pres^'nteu in the previous chapters. 

4*2 Instrumentatio n; 

Small scale vibratory tests usually employ mechanical 
oscillators which generate sinusoidal forcing function. Though 
the generation of sinusoidal forcing function is not an absolute- 
necessity, analyses of the data become extremely complicated 
With any other type of forcing functions. Again, any type 
of function, used in the analysis, may be transformed into 
an equivalent sinusoidal function. 

In the present in/esti_,ation an electrically 

operated oscillator, under the trade name “Unholtz-Dickie" * 
Shaker was used. 

* Vibratinj testing System No. 72 

Shaker: ^odel No. lo5, by Unholtz-Dickie Corporation. 

Capacity: Generated force: 180 lbs. peak 
Maximum free table acceleration; 42 g. 

Shaker Stroke: 1 inch 

Free table resonant frequency: 5000 cps . 
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The test equipment available here includes three 
separate vibration Test Systems. Each system is rated at 
180 lbs of generated force. Usable frequency range is 5 cps 
to 100 Tecs . Each system can be operated independently of the 
sytetems can be driven from a common oscillator. 

4.2.1 Simultaneous Operation of Three Systems: 

In the system of operation as shown in Fig. 4«1/ tli® 
console containing the BK 1039 Servo-oscillator is known as 
the master console. The t^^o consoles without oscillators are 
sTave consoles. In this mode of operation all Shakers are 
driven simultaneously at the same frequency. The BK 1039 
servo-oscillator maintains the acceleration or displacement 
level at one point constant. 

The arrangement/ shown in Fig. , will increase 
ofdecrease the vibration leve] of each shaker simultaneously 
by an equal amount, if the vibration level of each shaker 
varies* E^qh shaker is equiped with i .s own BK 1039 seivo-oscd.1 1- 
ator. In this case the 1039' s are interco mected in such a manner 
that the oscillator section of the 'Master' 1039 drives all 
the shakers, but the soivo-oscillator section in each 'slave' 

103 9 maintains the acceleration level at the corresponding 
shaker at the desired level. 

4.2.2. Single Operation: 

Each system can be operated independently. The slave 
console requires an oscillator signal input. 
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The operation employed here is, the master is directly 
attached to one of the three shakers and acceleration level 
is made conscant, with varying freqaency from 25 cps to 250 cps . 

4.2.3 Description of System: 

The vibration system consists of three parts, the 
shaker, the control console, and the Electronic Power Amplifier. 

4.2.3a Shaker: 

The shaker is tho important part of the system. It 
provid.-s th-- vibrati ng table to which the items to be tested are 
attached. The principal parts of the shaker is shown in the 

Fig. l\.2. 

The water cooled field coils energize the shaker body 
with flux paths as shown creating a high magnetic flux density 
in the annular air gap. The field coils are connected to a 
d--c power source. The armature assembly, consisting of the 
table and the driver coil, is suspended from the shaker body 
by the Unholt z-Dickie table suspension. Whan an alternating 
current is passed through the driver coil, an electro dynamic 
force is generated in the driver coil wires causing the table 
to move as indicated by the double ended arrow drawn on the 
table. The frequency of the alternating current determines the 
frequency of vibration in the table. The motion ofthe table 
is limited by the maximum stroke of the shaker, available 
voltage from tho electronic power amplifier m the low 
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frequeacics and by the rated force or maximum acceleration 
in the higher frequencies* The shaker is trunion mounted 
allowj.ng it to be oriented for vertical or horizontal vibration 
or any position between the limits. The operation employed 
here only thi- horizontal position of the shaker is utilized. 

4.2.3b Electronic Povjcr Amplifier: 

The electronic Pover Amplifier input is connected to 
the controlling oscillator signal source and the out put is 
connected to the driver coil of the shaker. The oscillator 
generates a low power, low voltage signal which is amplified 
by the Electronic Power Amplifier to the necessary power and 
voltage levels to cause the alternating cirrents to flow 
through the driver coil as required to generate rate force. 

The Amplifier is of straight forward design to provide the 
necessary cower levels through the operating frequency range 
without introduction of significant distortion. 

4.2.3c Control Console: 

The control console contains the Automatic Shaker 
Control, the Adjustable Power supply for the shaker field, the 
Dial-A-Gain and the control circuitry, pilot-lights, driver 
coil and field coil ammeters, interlocks and the system ' on-off ' 
switch . 

The Automatic shaker control consists of the 
Oscillator, the Auuomatic Level control and the Vibration 


Meter. 
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The oscillator originates the signal voltage which 
IS fed to the Power Amplifii. , which in turn causes an 
alternating current to pass thr' ugh the driver coil of the 
shaker . The driver coil current causes the shaker table 
motion. Thus the frequency and magn' tude of shaker table 
motion ss controlled by the frequency and voltage generated 
by the oscillacor, 

Tbe auton icic Jf-^l control adjusts the magnitude on 
che oscillator ou i put noltage fed to the power amplifier to 
account ■f"or van? lotis jn sn-tem frequency response due to 
resonant chara ?tP '"ictics of the items mounted on the shaker 
table and the inherent shaker sjstem response itself. 

In the ca' e here^ the freauency is varied manually 
while the automata c level conrrol maintains the set value 
of acceleration. 

The adjustable field suppily is provided to give more 
flexibilicy m system operation If fill rated generated 
force IS not required for a particular test, the operator 
may reduce the fi?ld cu’^rent in order to minimize the 
already lovj magnetic leaT'acc flux it the table surface or 
increase the obtainable lew frequency table displacement. 

Thfc Dial-^-Gain is a multipurpose instrument which 
allows monitanng acclcration from any piezo electric 
accelerometer. It includes a high impedance cathode 
follower input and a variable g-in amplifier. 
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4 . 3 Testing Procedure; 

For the present test program, an aluminium tank of 
29.8 by 29.8 ems in cross section and 49.6 ems high is filled 
with s and (Grain size analysis of the sand is shown in 
Fig. (5. 12)). The t-ank is mounced on the base plate (50 cm x 50 
ems) v;hich inturn is fixed rigidly to the shaker table 
(Fig. A '3) . This shaker table can be subjected to vertical 
Vibrations ol a selected amplitude and frequency. 

A mil'l steel fraae is fabricated from which a metallic 
(brass) sphere of 2 cm dia. is placed on the soil surface. 

A lo id IS imposed on this sphere by means of a thin brass 
(3 mm dia) rod. There is one counter weight which is suspended 
from a f riction-less pulley by a non extensible silk thread. 
The other end of the thread goes to the small platform, through 
another pulley, from which the load is imposed on the sphere. 
The counter weight slided over a smooth square section rod 
which is graduated. The combination of the rod and the 
counter weight makes a vernier arrangement, as result, it is 
possible to measure the displacement ofthe counter weight 
upto 0.001 ft. (The arrangerr'cnt is shown in Fig. A -3) . 

Before running the each t<-st, the sand is compacted by 
vibrating the soil. The frequency level is increased gradually 
keeping the acceleration level constant. The sinking of the 
sphere at different loading condition and frequency is noted 
by noting th'- displacement of the counter weight. 
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For each set of run the acceleration level is kept 
constant at l.Og and the frequency is varied from 25 cps to 
250 cps. 

Using Visicorder the vibration characteristics of 
sand 15 recorded and typical responses are shown in Figs 

The experiment carried out with the soils of different 
Grain sizes (Fig, 5,12) . 



CHAPTER - V 


TEST RESJLTS 


5 . 1 Test Data; 

The experiment was conducted at frequencies ranging 
from 25 cps to 250 cps . At each frequencies for different load- 
ing conditions from the upoer disc of the loading platform, the 
sinking of sphere is noted from the vernier scale arrangement. 

The loading conditions are ranging from 1 lbs to 21$ lbs from the 
top. The sinking of sphere versus time relations are shown in 
the figures 5.1 to 5.4 and then 5.S to 5.11. 


The displacement-time and acceleration-time responses 
are noted by Visicorder and are shown in figures 5.16 to 5. 19 
The experiment a s carried over on two soil samples 


(Fig. 5.12) . 

5 . 2 Soil Behavior; 

The soil medium chosen is medium fine dry sand. The 
physical properties of nhe soil are: 

Sample I : 


Specific gravity, 
Minimum voi^l ratio, e 

Bulkdensity, Vg 


2.50 

.468 

1.74 X 10^ kg/m^ 


Sample II: 

Gg = 2.55 
e = .412 

^ = 1.81 X 10 ^ kg/m^ 


109 pcf . 


113 pcf. 


Specific gravity. 
Minimum void ratio. 
Bulk density. 
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The Grain size distributions are shown in the Fig. 5 iO. 
To obtain minimuin void ratio the sand is vibnated for one 
minute at each frequency from 25 to 250 cps and after such an 
operation the void ratio is calculated (Before vibration the 
void ratios were, .724 and .546, respectively) . According 
to M .1 .1 . classification the sand for the first set is a mediv wi 
sand and the second set falls under coarse sand group. 

5 .3 Vibro-viscositv: 

The experiment shows that when the sand sample is 
subjected to intense vibration the forces of internal dry 
friction between the particles vanish and assume the 
mechanical properties of a viscous fluid. Thus the solid 
bodies placed on this aand will sink into it with a certain 
velocity if their densities exceed the density of the sand. 

This mechanical property is designated by the coefficient of 
vibro-viscosity . 

The mam object of rhe experiments carried out 
so far is to determine this coefficient and apply this 
information to the sinking behavior of piles due to liquef actio/i 
First of all the viscosity will be determined using Stokes 
classical method where by neglecting any effect of vibration 
of sand and neglecting the inertia and history integral of 
the Eq, (2.1) . And then the actual solution for an oscillatin 9 
fluid will be brought in to justify the first one. 
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5.3.1 Stokes Solution; 

If a spherical body descends through a semi- infinite 
VISCOUS fluid with constant velocity V, Stokes law says, 

the VISCOUS drag, STT^dV, is equal to the buoyant 
force (ms - rnf) g. The terms have already been explained 
in the Chapter III . 

Thus, Downward force = buoyancy + resisting force 
Than is, T +'^- . = — Pf g + BtryUd, 

3Tr^ dV = -IT d^ ( fg/ - ff) g + ^ ... (5.i) 

Now, in addition to this it will be assumed that 
th^ loll on the upp^ r disc will only increase the density 
of the sphere in the following way; 

Let^F be the load acting on the sphere in addition to 
its own weight, then iC equivalent increase of 

unit weight then^ 

>1 ' _ 6F__ ... (5.2) 

’ E " TT d3 

and actual ^ will be, 

s 

= yi+9„ ...(5.3) 

where, 'J i is the unit weight of sphere material. 

• s 
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Now, from equation (5.1) 



t - Vf) 
18 V 


... (5,4) 


Usiny equation (5.4) and utilizing the velocities tot 
different cases es noted in the Fig. 5.1 etcthe vibro— viscosit)/ 
versus frequency relations are obtained 'Figs .5 .12 & 5.13) . 

To obtain a representative value of the Kinematic 
Viscosity, , the follov7ing deductions were used i 

Obtained ^ in Ecj, (5.4) is in che unit gmf -sec/cm'^ 
To convert this into Ibf - sec/ft^; 

gmf - sec 2 

1 _ .205 sec - Ibf/ft ... (5.5) 

cm^ 

• Now to get density of the sand, if, Sp. gr. = G , 

s 

density = 1.94 Gg Slug/ft^ 

Thus Kinematic Viscosity of sand, using Eq. 5.4, is 

X .107 , ft^/sec. ••• (5.6) 

s 

The values of V with various frequencies are given 


in the Table 5.1. 
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5,3 ,2 Vcx"! f iC3.t ion of Exp© mitiG nt a.1 nssultsi 

To verify the rcsalts obcained in ths previous section 
with the solution proposed in the Chapter III, the following 
compute cions are made; 

The Eq. (3.1) can be rewritten in the form: 

X = P2/ ^^2 + . X ... (5.7) 


in which, 

'f ~ Kinematic viscosity of fluid, d = diameter of the sphete 
R = Df nsiLy ratio; k = /^dded mass coefficient 
P 
P- 


= lS/(R+k) , Q = 

3_= P. + Q/T'T' 


1 + k 


^2 " R + k 


A = w a,^ Cos wt 
1 ° 


O 

= w a^ Sin wt 


( (R-l) /(R+K) ) g ; g = acceleration due to gravity 


Notations ire similar as they were described in the 
Secti in 3.2. 

Ef|. 5.7 IS evaluated using the same numerical integratior 
technique, using, as variable for different set of results. 

7in a loroximato value of can be obtained from equation 5,6. 
Now, no the variables R, d, w, g, a^ are known before hand, the 
only v-^ri-ibic remains to be determined is V . Taking 
neighbouring v ilucs of ^ , as it is ob ained from Eq. 5.6, the 
variation ot the results are reported by varying the values 
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ot-j . The results are shown, along with the experimental 
data/ in tht Figs. 5.14 and 5.15. 

To illustrate the utility of the proposed solution 
^■^’1* ^ • l) / on 8©t of tost results are considered and compared 
with i:h.- solution of Eq. 5.7, with. 

Applied load, F = 2.5 lbs 
Frequency, f = 200 cps 

So. Gravity, = 2.5 , and with other 

varinbh ; r-.maLnin-i constant. These results are plotted in 
the Fifj. 5,14. It can be observed from the Fig. that the 
Kin‘.-matic viscosity,')?, as obtained from this analysis is much 
lowt'r th:in the om.* calculated from Stokes solution (Eq. 5.6) . 

For ex.’jmolo, from the Fig. 5.14, the Kinematic viscosity 
range at which the experimental data falls is from 20o to 300 
ft^/soc, although the Kinematic viscosity is equal to 420 ft^/sec 

using Stokes solution Bq. (5.6), which is the same as was used 

1 

by Barkan"^. 

The reason for this deviation is due to the fact that, 
in obco.ining d? from Eq. 5.6, the oscillation of the fluid is 
not tain.n into account and also the solution, is an approximate 
one unin<;; steady vdocity of a sphere in a stationary fluid. 

Thin ar*'!' !!■ c nt: is applicable, as well, for the other set of 
(:;xpCTii;< v-ilu* s. 

It can b’i noticed from the figures 5.12 and 5,13 that 
there exists a definite frequency at which Vibr-viscosity is 
minimum for a particulars soil. In the case I, it is 210 cps 





and for Case II its very near to 170 cps . 

The mam draw back of the experiment is that the wall 
effect of container is totally ignored, although it may 
contribute some discrepency to the results obtained with the 
theory . 

5,4 Vislcorder Records? 

* To evalua+e the respcuse m the soil medium a 

magnetic pick up a s fixed on a platform in the soil medium, 
lu IS observed in the F 2 .gf. 516 to 5.19, that the response 
gives suLiacient accurate displacement— time relationship at 
lower frequencies but as the frecaicncy increases the response 
nature bi-corrcs some what erratic. Although the input 
acceleration in all the ceses where g, which were mentioned 

i 

in the Figs, 5.1 to 5.4 and then in 5,8 to 5.11, but the 
visicor'ier records show some deviations, hs for example, 
in the Fig. 5.4, constant acceleration is shown as g but 
visicorder records (Fig. 5,19) show that the magnitude 
of acceleration is 0.7g instead of g. The calculations 
were made using the later acceleration, in comparing the 
laboratory results vjith the theory. 



TABLE 5.1 


Sample I : 


G =2.5, e 
s ' min 

= 0.468, 

Vg = 109 pcf 

Load 

(lbs) 

X 

X 

Frequency X 
(cps) 1 

Coefficient of | Kinematic 

viscosity,/!,, ^ viscosity,>?, _ 

( lbf-sec/ft2) xlO'^ ^ ( ft^/sec) x 10^ 

1.0 


25 

4.16 

0.857 

1.0 


30 

6.23 

1.280 

1.0 


50 

3 .12 

0.640 

1.0 


100 

2 .97 

0.610 

1.0 


170 

0.69 

0.140 

1.0 


200 

0.56 

0.117 

1 .0 


230 

0.62 

0.120 

1.0 


250 

1.25 

0.250 

1.5 


200 

1.31 

0.270 

1.5 


230 

1.31 

0.270 

1,5 


250 

0.51 

0.100 

? .0 


loo 

3 .04 

0.620 

2 .0 


200 

1.74 

0 .320 

2 .0 


220 

4.05 

0 .830 

2.0 


250 

3 .47 

0.710 

2.5 


190 

1.89 

0.390 

2.5 


200 

0.20 

0.042 

2.5 


220 

1.16 

0.230 

2.5 


250 

3.78 

0.770 
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Sample II : 

Gg - 2 . 55 , ©min 

= 0 . 412 ^ = 113 

pcf . 

1 

Load 1 

(lb) ' 

1 

_ 1 

T 

Frequency ] 

( cps) • 

f 

Coefficient of 

viscosity, 

(Ibf-sec/ft^) xl5‘^ 

' Kinematic 
' Viscosity, 'f , 

I ( ft^/sec) xio'^ 

1 .o 

170 

1.04 

0.21 

1 .0 

190 

1.04 

0.21 

1 .0 

200 

1.56 

0.32 

1.0 

250 

1.56 

0.32 

1 . 5 

170 

1.02 

0.21 

1.5 

190 

2.30 

0.47 

1.5 

200 

2.30 

0.47 

2.0 

170 

2.02 

0.41 

2 .0 

190 

2.02 

0.41 

2.0 

200 

2.43 

0.50 

2.0 

250 

4.05 

0.83 

2.5 

190 

1.89 

0.39 

2.5 

200 

2.04 

0.42 

2.5 

250 

2.04 

0.42 


CHifiPTER - VI 


GENERJiL DISCUSSION RND CONCLUSIONS 


Exp^-rimental and analytical studies have been 
conducted in this investigation to study the motion of a 
sphere in a vibrating soil medium. The vibro-viscosity of 
the fluidised soil medium has bean evalua^ied for a wide 
range of frequencies using Stokes solution and some of these 
values are compared with the modifit-d analysis presented in 
Chaoter III . The effect of gram size frequency/ load applie 
on the sphere etc. on the accelerated motion ofthe sphere 
have boon studied and are graphicilly presented. Relevant 
conclusions have besen presented, earlier, for each chapter. 
General conclusions from this investigation have been listed 
belows 

The experiments conducted, here in, reveal that, 
ev' n if the ground acceleration is lower than that of the 
value given by Barkan^ as threshold acceleration (l.5g) for 
fine to m-dium grained dry sand, liquefaction during the 
process of vibro floatation can occur at a lower value of 
acceleration ( .66g) if the frequency of vibration is high. 
This also depends on the duration of motion. This verifies 
D‘ Appolonia' s^*^ observation on vibro-floatation. From 
rhe experimental observations, made earlier, it can be 
observed that there is a particular frequency of the soil 
motion at which the resistance to the sphere penetration is 
the least. The cases observed here, for the first sample 
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thLs critical frequency is 210 cps an''! for the second 
case, it IS 170 cps, for the dv^nsity an-l void ratio 
mentioned there m. It is worth mentioning that at lower 
frequencies the top surface of the soil behaves like a 
susp-nsir,n upto a certain depth through which the ball can 
penetrate at a faster rate (Fig, 5.1) but fcheo at becomes mor' 
or less stand still (velocity being extremely small) . Thus 
th . criterian for liquefaction during vibro-f loatation can 
not be evaluat--l only emohasizing on critical acceleration 
but It shoul 1 includo rho frequency and the duration of 
vibration, as well. 

Based on the above informations, it is observed that 
a quantitative value of the soil behavior, namely kinematic 
viscosity of tht. soil, can be obtained. This kinerr.atic 
viscosity Will bo useful in predicting the behavior of 
structural elements founded on soil as a result ofsoil 
vibration. In obtaining a quantitative value of the 
viscosity terms, the classical m^ thod^ gives uncertain 
values which can not be correlated with the soil vibration, 
iis such it depends only on the sinking velocity and indep- 
endent of rhe nature of soil vibrarion. The proposed 
solution IS Jir ctly from th*- Navier-Stokes equations 
negltictinq convective acceleration terms and, hence, is 
more representative to the actual situation. Of course, 
tor t-h.^ presenc situation, where viscosity terms are much 
higher ir comparision to the inertial terms in the 
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Navier-Stokes equations, the acceleration, term is, a bit, 
insignificant but vibration ofthe soil medium has to be 
taken care off. 

Theoretical investigations, presented in these 
studies have the following advantages; 

1) The solution procedure reported, herein, for the 
accelerated motion of a sphere in an oscillating 
fluid, IS sim,o3e an 1 can be applied for computations 
for a wide ranae of parameters involved. 

2) These results can be used for practical solutions 
of sedimentation in reservoirs, grain size analysis 
end evaluation of the parameters involved such as 
viscosity, density ratio etc. 

Having established the stability of the numerical 
integration technique it can be applied for a general 
type of such problems. 


3 ) 
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SUBROUTINE EUNC. FIGURE 

HOTATIOFS ; 

N - 
H ; 

A,B,C 5 

1(1) t 

Y(2), Y(3): 

Q(1),Q(2) • 

DY(1),DY(2) • 
etc. 

R = Densit^r Ratio; RK = Added mnss coefficient; 

RA = Amplitude of vibration to characteristic length ratio. 

RB = g/vv^L = x'icceleration ratio 

PM = ((R-1 )/(R+RK)) * RB/n 4 

s 

N = SQRT(HEU/w*d**2) = Equivalent Stokes number 

3 

p = i8./(r+rk:) 


NO. BIPFERENTIdL EQUATION = 2 
STEP SIZE OP lUDUPEIIDEFT VARIABLE 
COE'R’FICIErx'3 IN T^ RUITGE-KUTT/i IffiTHOD 
IKDEFENBZl^T V'^RI^PLE 
DEFEI'BEIIT VARIABLE 
COERE'TIOj'i Ut-lIhI) TO EAiH V..RIABLE 
DIFPERET’TIAL EQUATIONS 












APPENDIX - II 


aocelerated motion of a sphere in an oscillating fluid 

2 3 

The original derivation is due to Basset ' , who 
formulated the problem for the accelerated motion of a sphere 
in a stationary fluid under the action of gravity. To suit 
the problem discussed earlier, the following computational 
effort IS made; 

» Det che sphere is of radius a, surrounded by 
a semi infinite fluid medium which is, let for the time being# 
at rest and sphere is moving with a velocity V. If squares 
and products of the velocity of the fluid are neglected, the 
current function must satisfy the differential equation: 

D(D — L~)y=o ...(1) 

dt 

where 

d^ Sin© d td \ ^ 

D + — (Cosecd — ) , and the (r, w ) 

(^3^2 r2 dO dB 

are polar co-ordinates of a point referred to the center of 
the sphere as origin. 

Now, in a single deg re e-of -freedom system, the velocity 
distribution function V can be written as, V = k - a, if the 
oscillation of the fluid be denoted as a = a^ Sin wt, the 
terms were similar, as described in the Chapter III. 




Let, R and 9 be the components of velocities of 
fluid along and perpendicular to the radius vector, then, if 
we assume that no slipping takes place at the surface of the 
sphere, the surface condition. 


R = 


Sin 9 




de 


= V Cos & 


9 = - 


a SmQ 


dvk 

d'ti 


- V Sin 6 , at infinity 

must vanish. 


If the impressed force is a constant force, here gravity, which 

acts in the dir ction of motion of sphere, and Z is the 

2,3 

resistance due to fluid ' , 


Z = 2 IT a 


It 


V 


(pa Sing- ^ Sin2 a) Sin© 


dS 


Can be written as. 


.TT 




U a 

^ dt 


+ m^g, 


(a + 2 \i; ) Sin^^dS 

dr ^ ^ 


O 


Vo' functions of 


r and t 


2,3 


m. 


Z = - 


^ { a +2 y, ^ 

dr ^ 


dt 
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The tertriy (a 


d y/. 
dr 


2.3 


+ 2 ''P 2^ reduces to ' 


- ^ +9a s/^t/TT 


+ ^ 

? ' 


Now, t has to be changed to C , V to F‘ (t-. 'J') dT , 
and integraeed the result with respect to T from O to t and 
obtained; with F(o) = O, where F(t) is the variable velocity, 
and F(o) = 0 gives rise to x = a, i.e. x 


Z =: 


0 


it ( f ‘ 

t 




t =o 


a2 


= a w, 
o 


TT + ^ + mfg 


m. 



( I + 9 ayT^ 


tr 


) F' (t -1> ) dr 
+ V +■ m^g. 

The equation of motion can be written as (Equating all the 
forces involved) ; 

,t 


** 1 t't 9 m 

X + tanr; + 

® f =2 


£ _1_ 

dt 


(u 


0 


(t- r ) ^C+ a ^ I dl^= mf 

+ (mg - mf) g 


... ( 2 ) 

where dots and slash denote derivatives with respect to t, 
and partial derivative with respect to t, respectively. 

Now, let us consider the integrals; 



66 


Integrating by parts. 


'■ (t -"t ) ) ^-5 + a \/'^^/-Tr =_p(t-T)(i 2 


't 


+ \ F(t -T) + -f- 


ihn + 

V%T 'I 


o 


d ]]' 


Putting F(0) = o, 
,t 


0 


Thus 


d 


J dt 


F{t -t;) a-?: (Sj-; + a- ) 


( 


t 


FCt-l) t + a- ^'P/n-i: j 


ai; 


0 


0 


'•(t-T) 


hf 


o 


t 

F'(t-T) dr + f \ ^9/jri^ F'(t~l) dT» 

'^o 


= V + ^ 

2 


^T'^/rr 


o 


F' ) 


s/ t -X 


dT 


Now, equation (2) bcco.iies , 
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Which 


ultimately reduces to; 
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The 


equation 4 .s the governing eguatton 
equatioQ 3.i of Chapter lii . 
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1 Reciprocoi of Vibro- Viscosity/ 
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VARIATION OF Ns. 
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Density ratios are shown in the bracket. 
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COMPARISON WITH SOME gXISTIMG RESULTS 
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INSTRUMENT INTERCONNECTION i 

(Reproduced from UNHOLT Z-DICKIE CORP Drowing No 9103269 
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